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0. Introduction 



One of the most fascinating applications of the Kac-Moody theory is the use of affine 

^ ■ Lie algebras and their groups to exhibit hidden symmetries of soliton equations. In 1981 

M. Sato [19] and Drinfeld-Sokolov [6] (see also [7]) discovered fundamental links between 

soliton equations and infinite-dimensional Lie groups. In an important sequence of papers 
"xT \ 

[2], [3], [4], Date, Jimbo, Kashiwara and Miwa gave a construction of the Kadomtsev- 

Petviashvili and Korteweg-de Vries hierarchies based on the representation theory of affine 

q^ • Kac-Moody algebras. The affine algebra a^ produces the KP hierarchy, while the KdV 



hierarchy is linked to s^. 

Date, Jimbo, Kashiwara and Miwa used the vertex operator realization of the basic 
highest weight module, which was discovered earlier by Lepowsky and Wilson [17]. In this 
realization for s^, the space of the basic module is identified with the polynomial algebra 
in infinitely many variables x = (x%, £3, £5, . . .). A completion of this space is interpreted 
as the space of r-functions. In this framework, the Backhand transformation that raises the 
soliton number of an N-soliton r-function, is given by the exponential of a vertex operator. 
The Casimir operator equation 

O(r®r) = (0.1) 

decomposes in a hierarchy of PDEs in Hirota form. 
The sine-Gordon equation 

u xt = sin(u) (0.2) 

is of great significance to both mathematics and physics. This equation appeared first in 
geometry and describes the surfaces of constant negative curvature in M. . The physical 
importance of the sine-Gordon equation is related to the fact that it is a soliton equation 
which is manifestly Lorentz-invariant. 
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The sine-Gordon equation played a special role in the development of the soliton 
theory. The Backlund transformation method, which is a precursor for the Lie theory 
approach, was first developed in the geometric setup for this equation. Before the con- 
nection between soliton equations and infinite-dimensional Lie algebras was discovered in 
1980's, Mandelstam wrote a paper [18] on quantized sine-Gordon equation, in which he 
used operators somewhat similar to the vertex operators. 

From the AKNS method [1] it is known that the sine-Gordon equation is related to 
sfo, but the representation-theoretic interpretation of this equation was missing. 

The main goal of the present paper is to fill this gap. Our starting point is the following 
observation. From the original Hirota's paper [13] we know that the sine-Gordon equation 
(0.2) can be written as a system of equations in Hirota form on a pair of functions (to, T\) 

D x D t {r Q oTQ -Tioti) = , Q ^ 

D x D t (T on) =t t 1 , 

where u is related to (to,ti) by u = 4arctanf jpj. To describe the soliton solutions in 
terms of r- functions, we to put To and T\ together into r = To + %T\. Then the iV-soliton 
solution for the sine-Gordon is given by (cf. [13]): 



N 




k=ll<j 1 <...<j k <N l<r<s<k V Jr Js/ \\s=l J \s=l 



Z js I X 



with real- valued a/s and Zj : s. 



If we compare this with the expression for the iV-soliton solution for the KdV hierarchy 

[11], [14] 




/ , Z js I Xm 



k=ll<J!<...<j k <N l<r<s<k •> J' J*< \ mF N„„ \ s = 1 



then it becomes clear that time t in the sine-Gordon equation is nothing but variable 
X-i missing in the KdV picture! The "correct" set of variables should include both pos- 
itive and negative odd integers: x = (. . . ,xs,X-i,xi,X3, . . .). We shall use the symbols 
^odd, ^ev,N dd, N ev to denote odd/even integers/natural numbers. 

Our extended hierarchy will be completely symmetric in positive and negative direc- 
tions. This means that the s/2-module that we need to construct should not be of the 
highest /lowest weight. However, the Casimir operator is a crucial ingredient for the con- 
struction of the hierarchy of differential equations (see (0.1) above). Formally, this operator 



is well-defined only for the highest /lowest weight modules. Nevertheless, we are still able 
to use the Casimir operator in our situation. 

In the vertex operator realization of the basic highest weight module, s/2 acts by 
differential operators on the algebra of polynomials C[xi, £3, X5, . . .]. In this paper we 
represent s/2 by the same differential operators, but now acting on the algebra of differential 
operators in variables xi,xs,X5, . . ., by left multiplication. The algebra of differential 
operators is isomorphic to the Weyl algebra W generated by Pi,qj, i,j G N dd with 
relations 

PiPj=PjPi, qiQj = QjQi, PiQj -QjPi- CiSij -1, Ci G C\{0}. 

Though the Weyl algebra is non-commutative, still the partial differentiations ^-, 
-Jf- are well-defined. This allows us to view the elements of the Weyl algebra as functions 
on a quantized space. 

The Weyl algebra has a Poincare-Birkhoff-Witt decomposition 

W = C[q 1 ,q 3 ,...]®C\p 1 ,p 3 ,..]. 

Since both subalgebras C[q\, qs, . . .] and C[pi,p3,...] are commutative, there is a linear 
bijection between W and the polynomial algebra C[. . . , xs,X-i,xi,xs, . . .]: 

TV : W — > C[xi,X3, . . .] <g> C[x_i,x_3, . . .]. 

The image of a differential operator under map tt is called the symbol of a differential 
operator (see e.g., [16]). 

The idea of our approach is to write a Kac- Moody group invariant equation on a Weyl 
algebra element f and then convert it into a hierarchy of partial differential equations in 
Hirota form on the symbol r = vr(f). Then we can use the Kac-Moody group action to 
generate the soliton solutions from the trivial solution f = 1. 

The function f = 1 is not anymore a solution of (0.1), since we have changed the 
action to be the left multiplication and 1 now represents the identity operator. However, 
f = 1 is trivially a solution of 

fl(f (8) f ) - (f <g) f)Q = 0. (0.4) 

This algebraic equation transforms into a hierarchy of the non-linear PDEs in the variables 
x = (. . . , xs,X-i,xi,xs, . . .) which contains two KdV subhierarchies going in positive and 
negative directions. 



The simplest new equation in this double KdV hierarchy is the following (cf. (4.34)): 

d 

Uyt = 7T- {U xxy + U x U y + U x ) (0.5) 

with x = x\,y = X-\,t = x%. This is a generalization of the KdV with two spatial 
variables, but different from the Kadomtsev-Petviashvili equation. 

The sine-Gordon equation can not possibly occur in the double KdV hierarchy (0.4) 
since the system of Hirota equations (0.3) can not be expressed in terms of r = r + %T\ 
alone. 

Using the boson-fermion correspondence we construct "skew" Casimir operators that 
split r into even and odd parts (Proposition 2.2). It remains a mystery whether the skew 
Casimir operators are related to the Dirac operator. We construct a hierarchy of equations 
by considering skew analogs of (0.4) (see (2.7), (2.8)). 

As expected, the sine-Gordon equation with respect to the variables x\ and X-\ ap- 
pears in it, "linking" the two KdV hierarchies. 

Here is a brief description of the structure of the paper. In Section 1 we construct a 
representation of affine Lie algebra a^ on a completion of the Weyl algebra W. We also 
review the boson-fermion correspondence and the construction of the Casimir operator. 
In Section 2 we consider the reduction of this representation to s/2, introduce the skew 
Casimir operators and study their properties. In Section 3 we show how to convert the skew 
Casimir operator equations into a generating series of PDEs in Hirota form and construct 
TV-soliton solutions for these equations. In the final section we derive non-linear partial 
differential equations from the Hirota equations. We also obtain the Leibnitz formula for 
the Hirota formalism and use it for the treatment of systems of Hirota equations. 

Acknowledgements. This work is supported by the Natural Sciences and Engineer- 
ing Research Council of Canada. I am grateful to Jack Gegenberg for his suggestion to 
study this problem. 



1. Boson-fermion correspondence and representations of a^ 

In the previous works on applications of infinite-dimensional Lie algebras to soliton 
equations the main idea is to represent Lie algebras by differential operators. In this setting 
the Casimir operator equation becomes an infinite hierarchy of non-linear PDEs in Hirota 
form. The fact that the Casimir operator commutes with the Lie algebra allows one to 
construct soliton solutions of this hierarchy. 



The approach that we take here is to use the action not on the space of functions, 
but rather on the Weyl algebra, or, equivalently, on the space of differential operators 
themselves. 

We consider the Weyl algebra W, an associative algebra with 1 generated by the 
elements pi, qi, i G N with the defining relations 

4 
PiQj - QjPi = — :%1> PiPj = PjPi, QiQj = QjQi, hj £ N. (1.1) 

The scaling factor —% in the first of these relations is chosen to put the vertex operators 
considered below, in a more symmetric form. Essentially, W is a factor of the universal 
enveloping algebra of the Heisenberg algebra in which the central element is identified with 
1. 

The Weyl algebra has a natural representation on the space of polynomials in infinitely 
many variables (Fock space) F = C[xi, x 2 , ■ ■ ■] by the operators of differentiation and 

multiplication: 

2 d 

Pj >-»• — r^— , Qj >-»• 2xj, j G N. 
J OXj 

Thus the Weyl algebra can be also interpreted as the algebra of differential operators on 
F. 

In spite of the non-commutativity of the Weyl algebra, the partial derivatives -J^— , -^- 
are well defined (which follows from the fact that the relations (1.1) survive the partial 
differentiation). These partial derivatives are inner derivations of W ([5], section 4.6): 

^- (/(q,p)) = \ fe,/(q,p)] , ^ (/(q,p)) = -\ fe,/(q,p)] , /(q,p) e W. 

(1.2) 
Here q= (q 1 ,q 2l ...) : p= (pi,P2,-.-)- 

The exponentials exp(a^-), a E C, (resp. exp(a^-)) are the shift operators 
Pi h ^ Pi + <%«, Qi l— ¥ Qi, (resp. pi i— >■ p^, ^ i— >■ qi + Sija), which are well-defined com- 
muting automorphisms of W. 

Clearly, an element /(q, p) G W can be written in many ways. We call the normal 
form of /(q, p) its presentation in which the generators ^'s are grouped to the left ofpj's in 
all monomials. The normal form can be found using the standard Poincare-Birkhoff-Witt 
procedure. It also corresponds to the decomposition W = W~ <8> W + , where W~ (resp. 
W + ) is an abelian algebra generated by qi, q 2 , . . . (resp. pi,P2, ■ ■ •)• 

The subalgebra W~ (resp. W + ) is isomorphic to the algebra of polynomials in in- 
finitely many variables C[x\, x 2 , ■ ■ ■] (resp. C[x_i, X- 2 , ■ ■ ■] )■ Consider the isomorphism of 



vector spaces 

n : W~ ®W + ^C[x 1 ,x 2 ,...X- 1 ,X-2,..], 

n (f(quQ2, ■ ■ -)g(pi,P2, •••)) = f(xi,X2, ■ ■ .)g{x- U x-2, . . .)• 

The map 7r coincides with the notion of a symbol of a differential operator expressed in the 
Weyl algebra language. This map allows us to convert the elements of the Weyl algebra 
into ordinary functions in commuting variables. 

We define the normally ordered product : fg : of /(q, p) and g(q, p) to be 

: /(q, p)#(q, p) := ^ _1 0(/(q 5 p)M#(q, p))) ■ 

We will consider a Z-grading of W — © W n by assigning degrees to the generators 

neZ 
as follows: 

deg(pi) = i, deg(qi) = -i, ie N. 

Next we introduce a completion W of the Weyl algebra, in which the vertex operators will 
be later defined. We set the completion of W n to be 

w n = n wr®w/" ( L3 ) 



i<0,j>0 
i+j=n 



and define W as 



W= © W n . (1.4) 

neZ 



It follows from the Poincare-Birkhoff-Witt argument that W has a well-defined structure 
of an associative algebra. 

In physics literature the generators of the Weyl algebra Pi,qi are called free bosons. 
The boson-fermion correspondence is a way of constructing free fermions out of free bosons 
and vice versa ([9], [14]). The free fermions are the generators t(ji,ip*, i G Z, of the Clifford 
algebra C^ satisfying the relations 

ipiip* + ip*ipi = 8 i:j , i/jiijjj + ipji^i = 0, ip*ip* + ip*ip* = 0, i,j e z. 

Form the formal generating series (fermion fields) 

i>(z) = J2 ^ z ~ k and !>*(*) = Yl ^ k . 



Proposition 1.1. (cf. [14], Theorem 14.10) The Clifford algebra C^ can be repre- 
sented on the space B = C[m, w" 1 ] ® W by vertex operators 



i/>(z) t-> z u d^uexp | - Y^ 1j zJ ex P 2 5Z p i 



z~ J 



ip*(z) i— ► w z " au exp , 



2^^ ex P _ 2^ Pj ' Z ' 



(1.5) 



Here the correspondence between two formal series is the correspondence between 
their respective components. 

Note that the components of the vertex operators belong to the algebra B = C[u, m -1 ]® 
W and act on B by left multiplication. The symbol z u 8^ is interpreted as z u 9^(u m ) = 



z m u m . 



The only difference between the statement above and Theorem 14.10 in [14] is that 
here we consider the action of these vertex operators by left multiplication on C[u, -u -1 ] <8> W 
and not on C[m,m _1 ] <E> F as differential operators. Proposition 1.1 is thus an immediate 
corollary. 

The obvious relations between the operator of multiplication by u and the fields 

^z),r(z) 

ip(z)u = zm/j(z), ip*(z)u = z~ 1 uijj*(z) 
can be rewritten for their components as 

i\) k u = uipk+u i>k u = u^k+i- (!-6) 

The classical matrix Lie algebras can be embedded into Clifford algebras and Weyl 
algebras. It is well-known that this can be also done for the affine Kac-Moody algebras [8] 
and, in particular, for the affine algebra a^ of infinite rank. 

The algebra cioo is a non-trivial one-dimensional central extension of the Lie algebra 
(Too of infinite matrices with finitely many non-zero diagonals: 



z2 a *3 E y 






3nGN \i — j\ > n =>- a,ij — 



GoO ^00 © (Li. 



The action of a^ on C[u, u 1 ] ® W is determined by the formula: 

i/jii/j* if i^j or z = j > 0, 



A;,^ ]p> 



" l -^i if i = 3 < 0. 



Using (1.5) we can write a generating series for this action ([14], 14.10.9): 



where 



»,ieZ Zl 



E ^^i^" J - y^h (7) "" r(, 1 ,, 2 ) - 1 



(1 °° ■ \ (1 °° 

2 E^l ~ Z 2>i J eX P I 2 E^ ~ Z 2~ J )Pj 

In order to get a connection with the differential equations in Hirota form, we should 
also consider the tensor square B ® B of the a^-module B — C[u, -u -1 ] ® W. 
The key link is the Casimir operator ([14], 14.11.2) 

keZ 



which acts on a completion B®B and commutes with the a^-action ([14], 14.11). To 
construct this completion we view W <E> W as a Weyl algebra in twice as many generators 
and take its completion as in (1.3), (1.4). The algebra B ® B is the tensor product of 
W <E> W with two copies of C[u, u~ x ]. 



2. Skew Casimir operators 

It is well-known that affine Kac-Moody algebra s/2 = s/2 (^[^t -1 ]) © CRT may be 
embedded into cioo by 

tf 2j+1 ^^£^ +2j+1 , Aj^^i-iy+lEij+j , ifnl , jeZ, (2.1) 
where 

The restriction of the representation of 0.00 on £> to this subalgebra yields (cf. [14], 14.13): 

k k 

H k ^~-p k , H_ k h^-q k , /ceNodd, (2.2) 



J]^-^I((-l)«*r(z)-l), (2.3) 

where 

r(*) = exp I £ ^ J exp I £ Pj z~i j . (2.4) 

Lie algebra g = sfo has a Z2-grading = 0o ©Si? where O is the principal Heisenberg 
subalgebra spanned by Hj , j E Z dd and K, while X is spanned by Aj + ^Sj^K, j e Z. 
The following lemma is an immediate consequence of the formulas (2.2)-(2.4) above: 

Lemma 2.1. (i) i« = ux for x G O , 

(ii) ot = —mx for x E Q 1 . 

This result can also be seen from the realization of elements of q as infinite matrices. 
The conjugation by u is a shift operator: u~ 1 E i jU = E i+ ij + i — SifiSj^l. It follows from 
(2.1) that the elements in q are represented by matrices invariant under this shift, while 
elements of g ± are represented by anti-invariant matrices. 

The subspace 1 <g> W in B is clearly invariant under a^ and s/2-actions. While the 
Casimir operator Q does not leave the space W (E)W invariant, it can be easily modified 



into an operator Aq, for which W <8> W is invariant: 

A = (l®«)0(ii _1 (8)1). 

Moreover using the conjugation by 1 (E) u k we obtain a whole family of skew Casimir 
operators: 

A fc = (l®M fc+1 )0(M- 1 ®«- fc ) , kZ. 

Note that the conjugation by elements u k <&l produces the same family because O commutes 
with u® u. 

It is easy to see that 

J2 ^z~ k = (ij(z)u- 1 ) ® (ur(z)). 

keZ 



Restricted to W® W, the generating series A(z) = Yl ^-kZ k can be explicitly writ- 

fceZ 
ten as 

A(z) = exp - Y^ <1j zJ ex P 2 5Z p J z ~ J ® exp ~2 ^ q i zJ exp ~2 ^ PjZ ~ 



The skew Casimir operators A k do not anymore commute with the action of a^, but 
as we see in the following Proposition, behave nicely relatively to the subalgebra SI2 in a^. 
Proposition 2.2. (a) For re G g 

(re (8)1 + l<8>re)A fe = A fc (re <8> 1 + 1 <8> re), kel. 

(b) For x G 0! 

(rr (8 1 + 1 <8) x)K k = -A k (x <8> 1 + 1 <8> re), he Z odd , 

(rr<8> 1 - l<8>re)A fe = -A fc (rc<g>l - l(8>re), keZ ev . 

Proof. Recall that A/- = (l(8>M fc+1 )0(ii _1 (8)M _fc ). Part (a) is obvious because xu = ux 
for re G 0o by Lemma 2.1 and re <S> 1 + 1 <8> x commutes with the Casimir operator Q. 

Verification of (b) is also straightforward. Suppose that re G g 1 and k is odd. By 
Lemma 2.1, re commutes with the even powers of u and anticommutes with the odd powers 
of u. Thus 

(re (8 1 + 1 (8 re)(l <8> m A:+1 )0(m" 1 <g> u~ k ) 

= (1 <g> M fc+1 )(rc <g> 1 + 1 <8> re)fi(u _1 <g> M~ fc ) 

= -(1 <g> M fc+1 )0(M _1 <g> u~ fc )(rc <g> 1 + 1 ® re). 

If k is even then 

(re <8> 1 - 1 (8 re)(l <8> M fc+1 )0(« _1 <g> u~ k ) 

= (1 <8> « fc+1 )(rc <8> 1 + 1 <8) x)0(m _1 <g> «" fc ) 

= -(1 <g> M fc+1 )0(M _1 <8 w~ fc )(rc <8> 1 - 1 <8> re). 

Let t = r(zi, Z2, • • •) be a formal Laurent series in finitely many formal real variables 
zi, Z2, ■ ■ ■ with coefficients in W. Let f*(zi, Z2, ■ ■ ■) be the complex conjugate of f. 
Proposition 2.3. The sets of solutions f(zi, Z2, . . .) of equations 

(f<g>f)A fc = A fc (f*<g>f*), keZ odd , (2.5) 

and 

(f <g>f*)A fc = A fc (f*<g>f), /cGZ ev , (2.6) 

are invariant under the transformation 

t(z!,Z2, ...)!->• f'(zi, z 2 , ■ • • , z) = r(z!,Z2, . • •) exp (aiT(z)) , 
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where agl and z is a formal real variable. 

The equations (2.5) and (2.6) are understood here as equalities of the corresponding 



W <® W-valued coefficients in the formal Laurent series. 

Proof. Since the components of T(z) represent the elements of g l7 we get from 
Proposition 2.2 that for k G Z dd 

exp (aiT(z)) <g> exp (aiT(z)) A k = exp (ai (T(z) (8 1 + 1 <g> T(z))) A k 

= A fc exp (-ai (T(z) 01 + 1® T(z))) = A k exp (-aiT(z)) <g> exp (-aiT(z)) , 
and in a similar way for k G Z ev : 

exp (aiT(z)) <g> exp (-aff(z)) A^ = A& exp (-aiY(z)) (8 exp («T(z)) . 

Suppose that f is a solution of the equation (2.5) 

(f<g>f)A fc = A fc (f*(8)f*), /cGZodd- 

Then for f' = f exp (aff(z)) we have 

(f ; ® f')A fc = (f ® f) (exp (aff(z)) ® exp (air(z))) A fc 

= (f <8 f )A fc (exp (-aff(z)) (8 exp (-air(z))) 
= A fc (f* <g> f*) (exp (-air(z)) <g> exp (-air(z))) = A fc ((f')* ® (r')*)- 

Thus f' is also a solution of (2.5). The case of (2.6) is completely analogous. 

Corollary 2.4. The series f = (1 + otiiT(zi)) . . . (1 + aNiF(zN)) in formal real vari- 
ables z\, . . . , zn with cti, ... , OiN G R is a solution for both (2.5) and (2.6). 

Proof. Trivially f = 1 is a solution for both equations. Thus by the above proposition, 
f = exp (a\iY(z\)) . . .exp (a^iT^N)) satisfies (2.5) and (2.6). Finally, exp (aiT(z)) = 
1 + aiV(z) because T 2 (z) = ([14], 14.11.15), and the claim of the Corollary follows. 

Consider the real and the imaginary parts of f : fo(zi, 22 > • • ■) = Re(f ), f\(zi, z%, ■ ■ ■) = 
Im(f). The equations (2.5) and (2.6) can be rewritten as the systems 



and 



(t <g> f - Tl <g> fi)A fc - A fc (f ® r - fi <8 h) = 
(f ® fi + fi fo)A fc + A fc (f <8 fi + fi (8) f ) = 

(f <8 f + fi (8) fi)A fc - A fc (f <8 r + fi <8 fi) = 
(to ® ti - fi (8 r )A fc + A fc (f <8 fi - fi <8 f ) = 
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fc G Z odd (2.7) 



fc G Z ev . (2.S 



The requirement that a and z are real may be dropped if we apropriately adjust the 
statement of Proposition 2.3: 

Proposition 2.5. The sets of solutions (fb,fi) of the systems (2.7) and (2.8) are 
invariant under the transformation 

(r ,r 1 )^(r ,r 1 ) = (ro,r 1 )^_ arW 1 j 

where z is a formal variable and a G C. 

The proof parallels the one of Proposition 2.3. 



3. Hirota bilinear equations and their solutions 

The differential nature of the Weyl algebra is encoded in its defining relations (1.1). 
In this section we shall see how to convert purely algebraic equations 

(f<g>f)Afc = A fc (f*<g>f*), keZ odd , 

(f<g>f*)A fc = A fc (f*<g>f), fceZev, 

into a hierarchy of differential equations in Hirota form on the "dequantized" functions 
To = 7r(f ),ri = 7T(fi). 

When a product of two "functions" in W is written in the normally ordered form, 
expressions involving partial derivatives (1.2) appear. One special case of this procedure 
is recorded in the following Lemma and will be later used in our calculations. 

Lemma 3.1. Let p, q be generators of the Weyl algebra satisfying pq — qp = c ■ 1. 
Then for f = f(q,p) we have 

and 

fexp(zg) = ex.p(zq) < exp I cz— ) r 

IV d Pj 

Here z is a formal variable and the above equalities are interpreted as equalities of formal 
power series. 

Proof. We shall prove the second identity, the first will then follow by applying the 
automorphism p i— > q,q i— ► —p. From (1.2) we get that 

f(q,p)q = qf(q,p) + c— = (? + c -^-J *(9>P)- 

12 



By induction we obtain 

f(q,p)q n = (q + c—j f(q,p). 

To complete the proof, we apply the above equality to the Taylor expansion of exp(zq): 

f(q, p) ex.p(zq) = exp I zq + cz— j r(q, p) 

= exp(zq) I exp ( cz— j f(q, p) \ . 

At this point we are ready to transform (2.7) and (2.8) into a hierarchy of Hirota 
bilinear equations. Note that the terms involved in these equations are of the form 
(/ <g> g)Ak ± A/-(/ (g) g). We shall work with such expressions using the generating series 
(f<8>g)A(z) ±A(z)(f®g), f ®g eW ® W. Denote the variables in the first copy of W by 
q'i,p'i and in the second copy of W by q" ,p'(- Then 

{f®g)A{z)±A{z){f®g) = 



oo \ / 1 oo 



M, P>(q", p") exp J - 2_^{q'j - q")z J exp I - 2_,{p 3 - p'-)z 

±ex P I 2^^' ~ q '^ zJ exp 2^^' ~ P 'J) Z ~' I /(*!'> p'Mq"* P'^ 
i=i / \ 7=i 

We shall assume that both /(q',p') and g(<i",p") are in the normal form. Moreover 
since p\ commutes with q'L we can replace f(q',p')g(ci",p") with the normally ordered 
product : /(q', p')g(q." , p") : , i.e. move all g's to the left of all p's. 

Making the change of variables pi = \{p'i + Pi), Pi = \{p'i~ p'DiQi = \Wi + Qi),Qi = 
\Wi ~ Q.'i)i we rewrite the above expression as 



/(q + q, P + p)#(q - q, P - P) : exp [ J^ q^z 3 exp ^pj 

k J'=i / \i =1 



z~ J 



±exp JZ^ exp ^PjZ J : /(q+ q, p + p)g(q - q, p - p) : . 

Collecting g's on the left and p's on the right using Lemma 3.1 (note that \fti,qj] 
4<$™ ■ 1), we get 



"1-3 
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oo 

exp ( ^2qjZ j I x 

exp | -2^y— J ±expf-2^^— — I I : /(q+ q,P + p)<?(q - q, P - p) : 

xexplj^z-' I . (3.1) 

We can convert this expression into the Hirota form. Recall that the Hirota bilinear 
differentiation is defined as 

P (D x , D y , . . .) [/(x, y, . . .) o ^(x, y, . . .)] := 

f d d \ 
P [dS' Q~T--)f( x + x 'y + y>-- )9{x -x,y-y,.. •) | i=0 ^=o,...- 

Using the Taylor formula we get that (see [14], 14.11.8) 

f d d \ 

p \^ z, d z, '"j x + x ' y + Vl ' ' ~ X} y ~ y ' ' ' = 

P(D x ,D y ,...) exp (xD x + yD y + . . .) [/(x, y,...)o g(x, y, ...)]. 
Applying this identity we transform (3.1) into 

exp [ J2 qj^ J i ex p I J2 v D <n 1 1 ex p ( _2 Yl j D Pi \ ± ex p ( _2 J2 z ~y Dq i 



exp ^pjDp. [: /(q,p)o#(q,p) :] I exp J^Pj^ J . (3.2) 

We use the reduction to the subalgebra sfo in 0.00 to construct solutions of (2.5) 
and (2.6). We see from (2.2)-(2.4) that this reduction involves only variables pj, qj with 
j G N dd, so the solutions we get in Corollary 2.4 do not depend on the even-indexed 
variables. Thus in (3.2) we may set D Pj = 0, D qj = for j e N ev - In (3.2) all g's are 
collected to the left of all p's and all g's are to the left of all p's. Because of that we can 
easily evaluate the image of this expression under the map it, which gives us: 

R±(z)[«(f)(x)oir(g)(x)], 
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where 



R±(z) = exp I 2_. x j z ^ ex P I /_. x jD 3 

VjeZ\{o} / \jeZ odd 



v-^ Z 3 » ' — ?~ J 



x ^ exp -2 J2 — D *-i ±exp -2 ^ —D Xj 

and x = (. . . , X-s, X-i, x±, xs, . . .). Here R±(z) is a Laurent series in z and a Taylor series 
in ajj. The coefficients of this series are Hirota differential operators. 

Now we can rewrite (2.7) and (2.8) as the generating series for the hierarchy of Hirota 
bilinear equations on tq = h(to) and T\ = 7r(fi): 

Res (z J R-(z)) (r oro-ri on) = 0, j 6 Z w , (3.3) 

Res (z?R+(z)) (r o n + n o r ) = 0, j G Z ev , (3.4) 

Res (^i?_(z)) (r o r + n o n) = 0, j'6 Z odd , (3.5) 

Res (V'.R + (z)) (r o n - n o r ) = 0, j E Z odd . (3.6) 

As usual, the residue denotes the coefficient at z~ x of a formal Laurent series. The coeffi- 
cient at each monomial in x's in the formal equations above is a Hirota bilinear equation 
on tq and T\. We shall study the Hirota equations that occur here in the next section. 
Now we turn to describing solutions of this hierarchy. 
From Corollary 2.4 it follows that 

t~o + ^Ti = 7T ((1 + aiiT(zi)) . . . (1 + a N iT(z N ))) 

is a solution of (3.3) - (3.6). From [14] we get that the right hand side converges to 

N , x 2 / / k 



i+y: e «*■■■«*** n (^^) H e (e y 

(3.7) 



fc=ll<ji<...<j fc <JV l<r<s<k ^ Zjr Zj ° J \meZ odd V s=1 



when \z\\ > \z2\ > ■ ■ ■ > \zn\- By analytic continuation we conclude that 
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[f] 

k=l l<ji<---<J2k<N l<r<s<2k 

2k 



Z j-r Z js 

Z jr "T Z js 



^ E EC ) X m\, (3-8) 



»eZ -- \ s=1 



v mfc£j odd 



[^] 



n= E ( _1 ) fe E a h--- a h k+ i n 

fc=0 l<Ji<---<j2fe + i<A r l<r<s<2k+l 



Z jr Z 3b 

Z jr "T" ^j s 



xex p( E Ec p™ ( 3 - 9 ) 

\m6Z odd V -=1 / / 

is a solution of the hierarchy (3.3)-(3.6) for all z\, . . . zn- 

As we shall see in the next section, this hierarchy contains the sine-Gordon equation, 
as well as two copies of the Korteweg - de Vries hierarchy. 



4. The KdV — sine-Gordon — KdV hierarchy 

In this section we transform the Hirota equations (3.3)-(3.6) into non-linear PDEs. We 
show that among other equations, this hierarchy contains the sine-Gordon equation, the 
KdV and the modified KdV equations. We begin by listing some of the Hirota equations 
that are coefficients at monomials in (3.3)- (3.6) (we multiply them by appropriate constants 
to avoid fractional coefficients). For simplicity of notations, we write Dj for D x . 

From the series (3.3): 

D-iDi{t g o t - ri o n) 

(D\ - D X D 3 ) (t o r - n o n) 

(D^Ds + 2D_ x Df - 3Dl) (to o r - n o n) 

(D-iAj + D\) (tq o t - n o n) 

From the series (3.4): 

(D_iDi - 1) (n o r ) 

(Df-DiD 3 )(rior ) 

(£>? - D_iD 3 ) (ri o to) 
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(at x±z), 


(4.1) 





(at xiz~ 3 ), 


(4.2) 





(at x_i^ -3 ), 


(4.3) 





(at x 3 z). 


(4.4) 





(at xiz), 


(4.5) 





(at xiz~ 3 ), 


(4.6) 





(at X32;), 


(4.7) 



(£>?-£>_!£>?) (no T ) = (atxlz). (U 

From the series (3.5): 

DI(t ot + n on) 
(^ + 2L» 1 D 3 )(rooro + r 1 or 1 ) 
(D_il>? - D_iD 3 ) (t o to + n o n) 
And finally from the series (3.6): 

(Df - D 3 ) (n o r ) 

(D^Df- Di)(noTo) 

(£>f - D 5 ) (n o T ) 

(2D* - 5D 2 D 3 + 3A.) (n o r ) 

One of the transformations that we consider below is u = 4 arctan ( — ) or, equiva- 
lently, ^ = tan (^). Let more generally ^ = ip(u), hence n = </?ro. 

In order to rewrite systems of Hirota equations as non-linear PDEs in function u, 
we need the Leibnitz rule Lemma for the Hirota differential operators. For a multi-index 
fi = (/3i,..., fi n ) let D 13 = D@* . . . Z)f™ be the Hirota differential operator and let d 13 = 
d^} . . . d@ n be the usual differential operator. For a pair of multi-indices fi and 7 with 

P\ (PA (Pn 






l at X-\Z j , 


(4.9) 





(at X1X3Z ), 


(4.10) 





(at x-\X\z~ 2 ). 


(4.11) 





(at x\z~ 2 ), 


(4.12) 





(at x-\z~ 2 ), 


(4.13) 





(at xiz -4 ), 


(4.14) 





(at x 3 z~ 2 ). 


(4.15) 



< 7, < Pi let 

1 7 J V 7i J V 7 
Lemma 4.1 



(a) Df*{<pfog) = E (^) drr^D^-^ifog). 

(b) D^/ o v^) = E (^) ^ 7 (^i ° ¥>2)D^(f o «,). 

The proof of this Lemma is straightforward. 
Sine-Gordon equation. Consider (4.1) and (4.5): 



£)_i.Di(tooto-Tioti) = 

(4.16) 
(D-iDi - 1) (n o to) = 

and set n = <^To- Applying Lemma 4.1, we transform this system into 

(1 - v? 2 )^-i^i(r o to) - D-iDifa o ^)r 2 = 
pD.iDifo o t ) + (a_i9i(^) - y?)r 2 = 0. 

Interpreting this as a system of linear equations on D_iDi(tq o to) and T5 2 we get that 

detf 1 "^ - D -^ l( r° P) )=0 (4.17) 
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since this system has a non-trivial solution. 

If we now substitute cp = tan (^) and multiply (4.17) by cos 4 (|-) then we obtain the 

sine-Gordon equation: 

f) 11 

sin(u). (4.18) 



dx-\dx\ 

If we let parameters otj in (3.8), (3.9) be purely imaginary then T\ will become purely 
imaginary, while To will remain real. In this case u will be purely imaginary, u = iU, 
■A- = tanh (-j), and the function U satisfies the sinh-Gordon equation: 



ITQ 

d 2 u 



wah(U). (4.19) 



dx-\dx\ 
Modified KdV equation. The same technique applied to (4.9) and (4.12) 

D\{tqo T q + Ti on) = 
(D\ - D 3 ) (n o r ) = 

yields 



(4.20) 



(4.21) 



(1 + <p 2 )Di(T o r ) + D 2 (<p o <p) T 2 = 

3d 1 (i P )D 2 1 (ro o r ) + (df(ip) - da{<p)){<p o v?)r 2 = 0. 

After the same substitution as above, cp = tan (l 1 ), multiplying the determinant of this 
system by cos 4 (j), we get the mKdV equation: 

du 1 / du \ d 3 



dxs 2 \dxij dx\ 



+ ^i (4-22) 



or for v = -jr*-: 

ax i 



dv 3 2 dv d 3 v 
dxs 2 dxi dx\ ' 

Mixed sine-Gordon - mKdV equation. The system of three equations (4.4), (4.7) and 
(4.9) 

{D\ + D_i£> 3 ) (to o r - n o n) = 

(^-D_ 1 D 3 )(r 1 or ) = (4.24) 

DKtoOTq + T! OTi) = 

becomes a system of three linear equations on D 2 (tq otq) , D-\Ds(tq oto), Tq and from 
its determinant 

I -if 2 l-tp 2 (D 2 + D-i£> 3 )(v? ^) 
det| <p -tp dlitp) - d^daiv) |=0 (4.25) 

1 + c^ 2 Df(<po(p) 
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after substitution ip = tan (f ) we obtain the equation 

d 2 u 1 . , . / du \ . , d 2 u , A rt „. 

9Z^=2 Sm(M) fej +COS(M) &f <4 ' 26) 

One can show that this equation follows from the sine-Gordon and mKdV equations (4.18), 
(4.22) taken together, however (4.26) may have an independent physical meaning. 

The second mKdV equation. The system of Hirota equations (4.9), (4.10), (4.14) and 
(4.15) 

(D\ - D 5 ) (ri o r ) = 

(2D\ - 5DlD 3 + 3D 5 ) (n o r ) = 

( 4 -27) 
D 1 (r O To + Ti ori) = 

k (£>? + 2£>i£> 3 ) (r o to + n o n) = 

applying the same method as above and eliminating x 3 using (4.23), is transformed into 
the following equation for v = 4 ~^- arctan ( ; 



1L 



dv d ( 3 r. 5 2 <9 2 i> 5 / dv \ d 4 v , , , . . 

^r^{~8 v + 2 v d4 + 2 v {d X - 1 ) + ^)' (4 - 28) 

which is essentially the second equation in the mKdV hierarchy constructed by the AKNS 
method [1]. 

The double KdV hierarchy. It can be easily seen that the hierarchy (3.3)-(3.6) is 
invariant under the transformation z i— > z~ ,Xj i— > £-j, i.e., it is symmetric in the positive 
and negative directions. 

If in (3.3) and (3.4) we set Xk = for k < and consider non- negative values of j, 
then 

Res (z j R±(z)) = ±Res (z j P(z)) , 

where 

P(z) = exp f^x jZ j exp ^ %£> x . exp -2 ^ — D x . . (4.29) 

V =1 / VieN odd / V jeN odd J / 

In this case the Hirota polynomials appearing in (3.3) and (3.4) differ only by sign. Com- 
bining (3.3) and (3.4) together we get that r = tq + %T\ satisfies 

Res (z j P{z)) (r o r) = for j e N ev U {0}. (4.30) 
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The Hirota equation 

(Dt-D 1 D 3 )(ro T )=0 (4.31) 



(cf. (4.2), (4.6)) after the substitution / = JU- ln(r) becomes the Korteweg - de Vries 



j_ 2 

equation 

9f , nt df , d 3 f 



8x3 dxi dx\ ' 

The series (4.30) is the KdV hierarchy [4], [14] in variables x±, x$, x$ . . .. By symmetry, the 
second copy of the KdV hierarchy extends in the negative direction x-i, xs, .... These 
two KdV hierarchies are linked together by the sine-Gordon equation (4.18) which involves 
variables x\ and x-\. The soliton solutions for the double KdV hierarchy are given by 
(3.7) with purely imaginary values of «i, . . . a at. 

Remark 4.2. Another way to construct this subhierarchy is to consider the Casimir 



operator fi a fj for affine Kac-Moody algebra sli (see [14]) acting on W®W. The double 
KdV hierarchy then arises from the equation 

(f<g>f)£2aff = Q aff (f <g>f). 

The double KdV hierarchy that we obtain here contains more than just two copies 
of the KdV hierarchy. It also contains equations that involve both positive and negative 
variables. The simplest such an equation on the function r = tq + ir\ is 

(£>_iZ> 3 - 4D_ 1 L>? + 3D 2 ) (to T )= 0, (4.33) 

which follows from (4.3), (4.4), (4.7) and (4.8). A similar equation (without the last 
term) appears in the D\ -hierarchy [15], but the solutions constructed here and in [15] are 
inequivalent. 

After the substitution g = -^— ln(r) we get 

9 2 g d ( A d 3 g nA dg dg n dg \ ,. n .s 

'4- ^ + 24— -2_— 2--3— 2- , (4.34) 



5x_i5a;3 dx\ \ dx-\dx\ dx-\dx\ dx\) 

which is a generalization of the KdV with two spatial variables, but different from the 
Kadomtsev-Petviashvili equation. 

It is interesting to note that the same equation appears in our hierarchy in a different 
way. If we combine (4.3), (4.4) with (4.9) and (4.11) then we get 

(LUD 3 - D_ X D\ + 3D 2 ) (r o r ) = (4.35) 
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and 

(LUDs - D- X Dl + 3DJ) (n o n) = 0. (4.36) 



These can be transformed into (4.33) by rescaling the variables Xj = 2 x 'j- Thus in addition 

1twt-i /"if I A e iA\ i\Tftr\ t- rniron r\\r I X / l lira ale 

9xi 



to the solutions g = ^- ln(r) of (4.34) with r given by (3.7), we also get solutions 



-=i+£(-d* e <•*-«*. n (ihrr) 2 

fc=l l<Ji<...<J 2 fc<A r l<r<s<2A; ^ Jr Js / 



X ex P , 9 



5 (E ^) *-i + 5 (E **.) *• + \ (E «?.) -3) (4.37) 



and 



T 






%r %'s 



c + 1 



fc=0 l<Ji<---<J2fe + i<A r l<r<s<2A:+l v Jr Js 



/2k+l \ /2fe+l \ /2fc+l \ \ 

x exp I - f J2 tf) ^-1 + 2 ( E ^ ) x i + 2 ( E z l) x A ■ ( 4 - 38 ) 

Miura transformation. It is well-known that solutions of the mKdV equation 

with V = jv (cf. (4.23) ) can be transformed into solutions of the KdV equation 

df df d 3 f 

dxs dxi dx\ 

by the Miura substitution 

The soliton solutions of (4.39) and (4.40) are given by V = 4^-arctanh ( i— J 

= 2g|- In ( ^^Yt ) an< ^ / = ^" m ( r o + i T i)i where To and T\ are given by (3.8) and 
(3.9) with purely imaginary parameters ctj. Note that tq ± %T\ is real in this case. 

A natural question to ask is whether for a given pair (tq, Ti), the functions V and / 
are linked by the Miura transform (4.41). The answer to this question is positive, and as 
we see from the following Lemma, the Miura transform is a feature of the whole hierarchy 
and not specific just to the pair of the KdV and mKdV equations. 
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Lemma 4.3. The functions V = 2 ^|- In ( ^°+^ 1 J and / = -^ ln(r + %t\ ) are related 
by the Miura transform (4.41) if and only if the Hirota equation (4.10) holds: 

D\{tqOTq + Ti ori) = 0. 

The proof of this Lemma is straightforward. 
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